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ABSTRACT 

We derive the first luminosity-metalhcity relation for a large objectively selected sample of local 
galaxy pairs and we compare the pairs LZ relation with the relation for the Nearby Field Galaxy Survey 
(NFGS). Galaxy pair members with small projected separations (s < 20 kpc/h) have systematically 
lower metallicities (~ 0.2 dex on average) than the field galaxies, or than more widely separated pairs 
at the same luminosity. There is a strong correlation between metallicity and central burst strength in 
the galaxy pairs. All five galaxies in the pairs sample with strong central bursts have close companions 
and metallicities lower than the comparable field galaxies. Our results provide strong observational 
evidence for a merger scenario where galaxy interactions cause gas flows towards the central regions, 
carrying less enriched gas from the outskirts of the galaxy into the central regions. The less enriched 
gas dilutes the pre-existing nuclear gas to produce a lower metallicity than would be obtained prior to 
the interaction. These gas flows trigger central bursts of star formation, causing strong central burst 
strengths, and possibly aiding the formation of blue bulges. We show that the timescale and central gas 
dilution required by this scenario are consistent with predictions from hydrodynamic merger models. 
Subject headings: galaxies:starburst-galaxies:abundances-galaxies:fundamental parameters- 
galaxies : interactions 



1. INTRODUCTION 

Galaxy interactions and mergers are fundamental to 
galaxy formation and evolution. The m ost wide l y sup- 
ported merger scenario is based on the iToomrd 11973) 
sequence in which two galaxies lose their mutual or- 
bital energy and angular momentum to tidal features 
and/or an extended dark halo and then coalesce into a 
single galaxy. Tidal interactions and associat ed shocks 
are thought to trigger star formation (e .g., lBushous3 
179871 IKennicutt et abllTflSTt iRa.rneFil2nn4l^. Theory pre- 
dicts that, as a merger progresses, the galaxy disks be- 
come disrupted by tidal effects, causing gas flows to- 
wards the central regions where kpc-scale starbursts may 
be fuel ed lB ajnesfcJIemaiiia3[12Qf3: iMihos & HernauTsfl 
IT9M llono et alJ 12004 ISpringelet al. 200^^ EvT 
dence for gas flows has been found in cold and warm 
gas in int eracting galaxies (e.g., ^ombcs et al. 1993 
Hibbard fc van Gorkom 1996: Georgakakis et al. 2003 



Chen et al.l,2002: Marziani et al...2003; .Rampazzo et al.l 
200,'Tt) . Large-scale gas flows combined with central 



bursts of star formation may be important for trans- 
forming merging galaxies into ell iptical galaxies (e.g., 
iKauffmann et al. 1993; Ko bavashil l2004t ISpringel et all 
ITOO,^- Nagashima et al. .200,C 
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The effect of galaxy interactjpns on star formation has 
been studied for decades. Mo rgan (.195^ first discov- 
ered "hotspots" in galaxies and IVorontsov- Ve lvaminovl 
(^959) published the first catalogue of interacting galax- 
ies. In 1966, Arp published the famous Atlas of Pe- 
culiar Galaxies, sparking mu ch research into th e prop- 
erties of peculiar galaxies. iSersic fc Pastorizal l|1967t ) 
found that the galaxies harboring hotspots all have ex- 
treme blue colors in their centers, a property of the 
majority of the galaxies in both Vorontsov-Velyaminov 
and Arp's catalogues. The suggestion that tidal forces 
in interacting galaxies could trigger bursts of st ar for - 
mation was first made by iLarson & TinslevI llT97l . 
and since then, numerous studies have provided ev- 
idence for inter action-induced starburst activity (see 
ISchweizeri I2?ffl5l for a review ) . Interacting galax- 
ies have elevated UV emission llPetrosian et all 119781) . 
Ha emission (e.g., 'Balzano"1983'; 'Kennicutt et al."198'J 
Bushousc 1987; Barton et al. 2000), IR emission (e.g., 
Lonsdale et al. 1984; Young et al. 1986; Heckman et all 
1986.: .gglgmgn fc Sage. 198^1. radio continuum emission 
(e^lHumme1in98lHCnndon et aUlTgSpjHummd et al. 
MW), and soft X -ray emission (e.g., iRead fc PonmanI 
1199 81 . Mergers have been linked to the extreme infrared 
luminosities seen in luminous infr ared galaxies (e.g., 
Soifer et aLlll984HArmus et alJlT987tlSanders fc Mirabel 



19961: iGenzel et al.lll998|) , and the unusually blue colors 
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observed in irregular galaxies ijSersic fc PastorizalllQSTfl . 

The observed relationship between star formation 
and galaxy interactions has been studied extensively, 
but we are just beginning to understand the metal- 
licity properties of merging galaxies. The metal- 
licity of normal disk galaxies is strongly correlated 
with galaxy mass. The first mass-metallicity rela- 
tion was found for irregular and blue compact galax- 
ies (jLeaueux et al.l 11979: Kinman & Davidson 1981). 
In subsequent work, luminosity was often used in- 
stead of mass beca use obtaining reliable mass esti- 
mates was difficult. iR.ubin et al.1 l) 19841) provided the 
first evidence that metallicity is correlated with lu- 
minosity in disk galaxies. Further work using larger 
samples of ne arby disk galaxies confirmed this result 
iBothun et al.l T984: Wvse fc S ilk 1985; Skillma n et all 
[1989: Vila-Cos tas fc Edmundsll9 92: Zaritskv et al .ll994l: 
iGarnett ■2002^') . The origin of the mass-metallicity re- 
lation for disk galaxies is unclear. A mass-metallicity 
relation will arise if low mass galaxies have larger 
gas fractions than higher mass galaxies, as is ob- 
served in local galaxies j McGaugh fc dc Blok 1997; 
iBell fc de Jcmel 1200(1: iBoselli et al...2001.) . Selective loss 
of heavy elements from galaxies in supern ova-driven out- 
flows can also account for the relation ijGarnettl 120021: 
l^emonti et al. 2004). 

In addition to providing information about previous 
episodes of star formation and mass-loss, metallicity may 
provide clues into the dynamical effects of interactions 
on gas flows in galaxy pairs. Isolated disk galaxies often 
displ ay metallicity gradients (sec Shields 1990; Dincrstcin 
IT9M for a review) which may be altered during a merger 
as a result of tidally-driven gas flows. Metallicity gradi- 
ents in interacting galaxies could potentially provide a 
unique method for finding interacting galaxies that have 
experienced strong tidal gas flows. Few studies to date 
investigate the effect of galaxy mergers and interactions 
on the metallicity properties of galaxies. Recent work 
in this area foc uses on small numbers of g alaxies, or on 
single galaxies. iDonzelli fc P astoriza' (200(f) analysed the 
optical spectra of 25 merging galaxies and find that the 
spectra of mergers are more excited on average than in- 
teracting pairs. Donzelli & Pastoriza attributed this re- 
sult to lower gas metallicity in the galaxy mergers, al- 
though it is not clear from their work what the cause of 
the lower gas metallicity could be. Maraucz ct al. (2002) 
used the [N II] /Ha ratio to study the metallicities in disk 
H II regions in 13 non-disrupted galaxies in pairs and a 
large sample of isolated galaxies. They find that the 
[N II] /Ha ratios in the disk H II regions of non-disrupted 
galaxy pairs are higher than for isolated galaxies, but 
that there is no difference in nuclear [N II] /Ha ratios . 
Recently, .Baldi et al. 1.2004) and Fab biano et all (j2004D 
investigated the Ne, Mg, Si, and Fe metallicities within 
the hot gas in X-ray luminous star forming regions of the 
Antennae galaxies. They find that the metal ratios for 
these hot regions are consistent with enrichment of the 
gas by Type II supernovae. 

These investigations provide the first clues into the 
impact of interactions on the metallicity gradients in 
galaxies. A more comprehensive analysis of the metallic- 
ity properties in merging galaxies requires well-defined 
and carefully selected surveys. In this paper, we com- 
pare the metallicity and star formation properties of two 



samples that have been objectively selected from the 
CfA redshift catalog: (1) t he galaxy pairs sample of 
iBarton Gillespie et alJ i2 00^. and (2) the Nearby Field 
Galaxy Survey (|jansen et al. 2000bJ). The objective sam- 
ple selection enables powerful statistical comparisons to 
be made between the two samples. We describe and com- 
pare the two samples in Section|21 In Section|21we discuss 
our classification scheme for the removal of AGN and we 
compare the locus of the galaxy pairs with that of the 
NFGS on the standard optical diagnostic diagrams. We 
investigate the ionization parameter differences between 
the pairs and the NFGS samples in Sectional The metal- 
licity estimates for the galaxy pairs and the NFGS are in 
Section [SI which contains the first luminosity- metallicity 
(LZ) for galaxy pairs. We show that galaxy pairs have 
lower metallicities on average than field galaxies within 
the same luminosity range. We investigate the relation- 
ship between metallicity and central burst strength and 
we examine whether metallicity is correlated with the 
presence of blue bulges for the galaxy pairs. In Sectional 
we discuss our results in terms of the merger scenario 
predicted by hydrodynamic simulations, and our conclu- 
sions are in Section |3 

Throughout this paper, we adopt the fiat A-dominated 
cosmology as measured by the WMAP experiment {h = 
0.72, = 0.29; Spergel et al. (2003)). 

2. SAMPLE SELECTION 

We use two ob jectively s elected local samples: the 
galaxy pairs of Bartoi^£t_^lJ_(^bOO) and the Nearby Field 
Galaxy Survey (Jansen et al....20 00b). Ground-based nu- 
clear optical spectra are available for both surveys; the 
spectra include ~ 10% of the total B-band light. 

Jansen ct al. (2000b) selected the NFGS from the CfAl 
redshift survey ( Davis fc PeeblesI [1983; Hu chra et alJ 
119831) . The galaxies were selected from 1 mag- 
wide bins of Mzwicky^ and CfAl morphological type 
to approximate the local galaxy luminosity function 
(e.g., Marzkc ct al. 119941) . To avoid a strict diame- 
ter limit, which might introduce a bias against the in- 
clusion of low surface brightness galaxies in the sam- 
ple, Jansen et al. imposed a radial velocity limit, 
VLG(kms-i) > io-" i9-0-2Mzwicky (.^^ith respect to the 
Local Group standard of rest). Galaxies in the direc- 
tion of the Virgo Cluster were excluded to avoid favoring 
a cluster population. The final 198-galaxy sample con- 
tains the full range in Hubble type and absolute mag- 
nitude present in the CfAl galaxy survey and is a fair 
r epresentation o f the lo cal galaxy population. 

iJansen et alJ l)2000al) obtained nuclear spectra with a 
3" slit centered on the nucleus and aligned along the ma- 
jor axis of each galaxy, sampling 3" x 6'.'84. The covering 
fraction of the nuclear spectra depends on the radial light 
profile and ranges between 0.4-72% of the light within the 
-626 isophote, with an average covering fraction of 10%. 
However for the sample that we use in this study (de- 
scribed in Section EJ, the covering fraction ranges from 
0.4 -33% with a mean covering fraction of 9%. 

We calibrated the nuclear fluxes to absolute fluxes by 
careful comparison with B-band surface photometry (de- 
scribed in Kgwl cv et al. 2002). We corrected the Ha and 
up emission-line fluxes for underlying stellar absorption 

^ Mzwicky is approximately equal to Mg 
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Fig. 1. — The B-band luminosity distribution of (a) the galaxy 
pairs and the NFGS, (b) the galaxy pairs and the NFGS with a 
recessional velocity limit cz > 2300 km/s, (c) the galaxy pairs 
and the NFGS that have measurable metallicities and are domi- 
nated by star-formation. The NFGS contains a larger fraction of 
lower luminosity galaxies than the pairs sample. For the purposes 
of metallicity comparisons, we compare galaxies within the same 
luminosity range of —21.5 < Mb < —18. 



as described in lKewlev et alJ l|2002|) . 

We used two methods to correct the NFGS emission- 
hne fluxes for Galactic extinction , base d on: (1) 
the HI maps of IBurstein fc HeilesI lll984|) . as hsted 
in the Third Refer ence Catalogue of Bright Galaxies 
l)de Vaucouleurs et al. 1991), and (2) the COBE and 
IRAS maps (plus the Lcidcn-Dwingeloo maps of HI emis- 
sion) of ISchlcgcl ct al. (1998). The average Galactic 
extinction is E(B - V)=O.OU ± 0.003 (method 1) or 
E (B - 1^)=0.01 6 ± 0.003 (method 2). 

iBarton et aP l|200nl hereafter BGKOO) selected a sam- 
ple of 786 galaxy pairs and N-tuples from the CfA2 red- 
shift survey, (|Gcllcr & Huchra 1989; Huchra etal]IT990l 
11995(1 . The CfA2 redshift survey is magnitude-limited at 
"mzw = 15.5. These galaxies were selected with line-of- 
sight velocity separations AV < 1000 km/s, projected 



separation s < 50 h^^ kpc, and cz > 2300 km/s. The 
lower velocity limit excludes the Virgo cluster and limits 
the angular sizes of the galaxies. Of these 786 galaxies, 
BKGOO obtained nuclear spectra for 502 galaxies using a 
3" slit centered on the nucleus of each galaxy. The cov- 
ering fraction of the spectra is approximately 10% of the 
galaxy light, comparable to the NFGS covering fraction. 
The spectra were flux calibrated in a relative sense by 
BKGOO using standard stars. The observations were not 
taken under photometric conditions, thus we have only 
relative calibration across each spectrum, suitable for line 
ratio diagnostics and equivalent width measurements. 

Because the BGKOO sample selects only on projected 
separation and line of sight velocity separation in a com- 
plete magnitude limited sample, the selection method 
is unlikely to produce correlations between pair separa- 
tion and any spectroscopic or photometric property. The 
BGKOO sample includes all pairs which satisfy the selec- 
tion criteria, regardless of morphology or gas fraction. 

The 502-galaxy sample is primarily composed of pairs; 
382/502 (74%) galaxies are in pairs, 70/502 (14%) galax- 
ies are in groups of three, and 10/502 galaxies (10%) are 
in groups of four or more. A total of 190 galaxies in 
the pairs sample have B-band and R-band magnitudes 
measured by BKGOO. 

Figure shows the B-band magnitude (Mb) distribu- 
tion of the pairs and the NFGS. The NFGS spans a much 
larger range of Mb , primarily because the NFGS sample 
was selected to reflect the local galaxy luminosity func- 
tion. Because of this selection, low luminosity galaxies 
appear to be over-represented in the NFGS sample com- 
pared to the pairs. The galaxy pairs were selected with a 
lower recessional velocity limit of cz > 2300 km/s. Fig- 
ure^ gives the luminosity distribution after applying a 
lower cz limit of cz > 2300 km/s to the NFGS. FigureUt 
shows that for galaxies with cz > 2300 km/s, the NFGS 
has a luminosity distribution that is slightly shifted to- 
wards high l uminos ities. This effect occurs because 
iJansen et al.l l|2000bD applied a lower radial velocity limit 
that increases with luminosity for the NFGS. For the pur- 
poses of metallicity comparisons, we compare galaxies 
within the same luminosity range of —21.5 < Mb < —18 
and with cz > 2300 km/s. This selection ensures that 
(1) we compare similar size galaxie s to minimise aper- 
ture effects (see lKewlev et al.l 120051 for a discussion of 
aperture effects in the NFGS sample), and (2) there 
are sufficient numbers of galaxies in each 1-Mag lumi- 
nosity bin for comparison. Because galaxy metallicities 
correl ate with global luminosity ( e.g., Kobulnicky et all 
199^; 'Melbourne fc Salzer' "200 1 iTVemonti et alJ 120041 
Kobulnickv fc Kcwlcv 2004; Sal zer et all 120051) . metal- 
licities must be compared within the same luminosity 
range. We use B-band for our sample comparisons be- 
cause both the pairs and the NFGS samples were selected 
from th e magnitude limited CfA redshift catalogs based 
on the jZwickv et alj|]'96ll) catalog. Therefore, the B- 
band is the consistent band to use for comparisons be- 
tween samples. 

3. CLASSIFICATION 

We use the nuclear optical emission-line ratios to clas- 
sify the galaxies into H II region-like or AGN-dominatcd 
classes. For galaxies with [O III]/H/3, [N II] /Ha, 
[S II] /Ha and (if available) [O I] /Ha ratios, we used 
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th e theoretical optical classification scheme developed 
bv lKewlev et alJ l|2001ar . In this scheme, galaxies that 
lie above the theoretical "extreme starburst line" in the 
standard optical diagnostic diagrams are AGN; those 
below the line are H II region-like. Galaxies in the 
AGN region in one diagnostic diagram but in the H II 
region section of the other diagram are "ambiguous" 
(AMB) galaxies. A fraction of galaxies in the NFGS 
(35/198) and pairs sample (49/502) have [N II]/Ha ra- 
tios but immeasurable [O III] or H/3 fluxes in their nu- 
clear spectra. We classify these galaxies as HII region- 
like if log([N II]/Hq;)< — 0.3, typical of starburst galaxies 
and H II regions re.g.. lKewlev et al.ll200ThL iDooita eraP 
120001 We are unable to classify those galaxies with 
-0.3 <log([N II]/Ha)< 0.0 if [O III]/H/3 is not available 
because these line ratios can be produced by both AGN 
and starburst galaxies. Galaxies without [O III]/H/3 but 
with strong log([N II]/Ha)> 0.0 are c lassed as AGN (e.g., 
Figure 1 of Brinch mann et a^l2004^ . The nuclear spec- 
tra of 124 NFGS galaxies contain sufficient line ratios for 
classification. Of these, the nuclear spectra of 105/124 
(85%) NFGS galaxies are dominated by star formation, 
12/124 (10%) are dominated by AGN, and 7/124 (5%) 
are "ambiguous" . Because ambiguous galaxies are likely 
to contain both starburst and AGN activity (see e.g., 
iKewlev et al. 2001b; HiU et al. 1999), we do not include 
them in the following analysis. The remaining galaxies 
do not have sufficient emission-lines in their spectra to 
allow classification. Many of these unclassified galaxies 
are ellipticals. Out of the 105 star-forming galaxies in 
the NFGS, a total of 43 galaxies have cz > 2300 km/s 
and —21.5 < Mb < — 18. These 43 galaxies comprise our 
comparison field galaxy sample. 

A total of 212/502 of the pair members spectra con- 
tain emission lines suitable for classification. Of these, 
153/212 (72%) pair members are dominated by star 
formation, 36/212 (17%) are dominated by AGN, and 
23/212 (11%) have ambiguous line ratios. Out of the 
153 star-forming galaxy p airs, a total of 92 galax i es hav e 
B-band photometry from iBarton Gillespie et al.l l)2003D , 
and 86 galaxies have -21.5 < Mb < -18. These 86 
galaxies comprise our "galaxy pairs" sample. Although 
all of the galaxies within our pairs sample are members 
of a pair or n-tuple, our final 86-galaxy sample does not 
necessarily contain all members of each pair or n-tuple. 
We note that classification (and metallicity determina- 
tion) with optical emission-lines selects against elliptical 
galaxies and galaxies with low gas fractions in the origi- 
nal 502-galaxy catalog of BKGOO. 

We check that our emission-line selected pairs sample 
has not biased the pairs against luminous high-mass (and 
high metallicity) systems as a function of projected sepa- 
ration in Figure 121 Clearly there is no deficit of luminous 
pair members at low projected separations. 

In Figure 0| we compare the position of the galaxy 
pair members with the field galaxies on the standard 
optical diagnostic diagrams. These diagrams are based 
on strong optical line ratios, and are co mmonly used to 
separ ate AGN from star-forr ning galaxies llBal dwin et aP 
IT98llVeilleux fc OsterbrockllT987t iKewlev et al.ll20oTbr 
Galaxies dominated by star-formation l ie below the 
theoret ical classification line derived by IKewlev et al.l 
l)2001a(l . The position of the star forming galaxies on 
these diagrams is determined primarily by (1) metal- 



T 



-22 




OC • • • 

-20 • • • 

-19 - - 

: : ?^ \ : : : \ : , : \ , : ^ 

20 40 60 80 

s 

Fig. 2. — Projected separation s in units of kpc/h versus R-band 
luminosity for the entire pairs sample. Filled circles denote the 
galaxy pairs that satisfy our emission-line selection criteria. Clearly 
our selection does not bias the sample against massive galaxies at 
low projected separations. 



licity, (2) the ioni zation state of the g as, and (3) star- 
formation history (IKew lev et alJl2001a|) . Figure |31 shows 
that the star-forming galaxy pairs l ie closer to the classi- 
ficatio n line than the field galaxies. iDonzelli fc Pastorizal 
l)2000|) found a similar result when they compared merg- 
ing with interacting galaxies on the diagnostic diagrams. 
They found that star-forming mergers occupy a higher 
locus on the diagnostic diagrams than well-separated, 
interacting galaxies. Donzelli et al. interpret this re- 
sult as a consequence of a greater level of excitation in 
the nuclear spectra of galaxy mergers and they suggest 
that the greater level of excitation could result from a 
lower gas metallicity in the mergers. High excitation 
is correlated with low gas metallicity; the fewer met- 
als in a galaxy, the fewer coolants are available and the 
gas maintains a high temperature. The high tempera- 
ture results in reduced recombination rates, and therefore 
higher excitatio n lines like [O III] are produced. However, 
IMiirauez et"an (2002) find that the nuclear [N II] /Ha ra- 
tios for 13 non-disrupted galaxy pairs are similar to those 
of isolated galaxies. We emphasize that the position on 
the diagnostic diagrams is not a simple function of metal- 
licity, and that the pairs locus may result from one or 
a combination of a lower metallicity, higher ionization 
state, or a lo ng (continuous) star formation history last- 
ing > 6 Myr ijKewlev et aU feOOla). In Sections H and [3 
we directly compare the ionization, metallicity and star 
formation properties of the galaxy pairs with the field 
galaxy sample. 

4. IONIZATION PARAMETER 

The ionization state of the gas is commonly described 
in terms of the local ionization parameter (q) , defined as 
the number of hydrogen ionizing photons passing through 
a unit area per unit density of hydrogen: 

9= (1) 

n 
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Fig. 3. — The standard optical diagnostic diagrams siiowing the 
position of the nuclear line ratios of the galaxy pairs (filled cir- 
cles) and the NFGS galaxies (un filled circles) w ith cz > 2300 km/s 
and -21.5 < Mb < -18. The IKewlev et al.l I2001alf theoretical 
classification line is shown in red. The two dashed lines indicate 
the error range (~ O.ldex) of the classification line (solid line). 
The star-forming pairs lie closer to the classification line than the 
star-forming field galaxies. 



where S^io is the ionizing photon flux through a unit area, 
and n is the local number density of hydrogen atoms. 
The ionization parameter q has units of velocity (cm/s) 
and can be physically interpreted as the maximum veloc- 
ity of an ionization front that can be driven by the local 
radiation field. Dividing by the speed of light gives the 
dimensionless local ionization parameter U = q/ c. 

Many of the common optical metallicity diagnostics 
are sensitive to the ionization parameter. So long as 
the EUV spectrum of the exciting source is reasonably 
well constrained, the ionization parameter can be deter- 
mined using the ratios of emission-lines of different ion- 
ization stages of the same element. The most commonly 
used ionization parameter diagnostic is based on the ra- 
tio of the [ O ml A5007 and [O III A3727 emission line 
fluxes (e.g-.. IKewlev fc Donita ll2002il) . Figure ^ shows the 
[O III]/[0 II] vs [ N II] /[O II] l ine diagnostic diagram and 
the model grids of' Kewlev etal.1 ll|20cilal|l . The mctaUicity 
Z is deflned in terms of log(0/H)-|-1 2 where the currently 
accep ted value of solar is 8.69±0.05 ijAllende Prieto et al.l 
120011) . Kewley et al. calculated these grids using a 
combina tion of the Pegase stellar popu lation synthe- 
sis code ijFioc R.occa-Volm crangc 1997) and the Map- 
pings photoionizat ion code (Sutherland & Dooita 1993; 
iGroves et al.ll20d^ . The models show that the position 
of galaxies on this diagram is determined by their ioniza- 
tion parameter, metallicity and star formation history. 
We examine two star formation history models: (a) a 
continuous star formation model, and (b) an instanta- 
neous zero age burst model. The galaxy pairs are shifted 
towards a greater [O III]/[0 II] ratio in Figure ^ than 
the fleld galaxies. Their position on these diagrams may 
result from either a greater ionization parameter, a con- 
tinuous star formation history, or a combination of the 
two. A greater ionization parameter could result from 
a greater number of ionizing photons per unit area or a 
lower hydrogen density. More ionizing photons may re- 
sult from a larger fraction of young, hot stars within the 
nuclear regions in the galaxy pairs compared with the 
fleld galaxies. This picture is consistent with numerous 
studies indicating that nuclear star for mation increases 
in galaxy interactions and mergers (e.g.. lBushouse et al.l 
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Fig. 4.— The [O III] A5007/[O II] A3727 ratio versus the 
[N II] A6584/[0 II] A3727 ratio diagnostic diagram. The galaxy 
pairs sample and the NFGS sample (with cz > 2300 km/s) are 
plotted as filled and unfilled circles respectively. T he grids show 
the locus of the theoretical star-forming models of IKewlev et all 
12)013) for two extreme burst scenarios: (left) continuous star 
formation, and (right) instantaneous zero-age burst. The metal- 
licity Z is defined in terms of log(0/H) + 12 where the currently 
accep ted value of solar is 8.69+/-0.05 dex lAUende Prieto et alj 
120011) . The ionization parameter q is in units of cm/s. The galaxy 
pairs sample contains more galaxies at relatively low metallicities 
than the NFGS. The galaxy pairs have higher [O III]/[0 II] ra- 
tios on average than the field galaxies, indicating either a higher 
moan ionization parameter, a longer episode of star formation, or 
a combination of the two. 



11988' 'Donzelh fc PastorizaM iggl IBarton Gillespie et alJ 
12003: Nikolic et al. 2004, and references therein). We 
note that shock excitation cannot explain the line ra- 
tios observed in Figures 0] and |3| Shock excitation raises 
the [O III]/[0 II] ratio, but it also raises the fS II] /Ha ra- 
tio by a similar amount fe.g. JSuther land fc Dopitalll993t 
IKewlev ct al. 2001b); we do not observe greater [S II]/Ha 
ratios for the galaxy pairs. 

We calculate the ioni zation parameter u sing the 
[O III]/[0 II] diagnostic in iKobulnickv fc Kewl ev (200^ 
using the metallicity estimates described in Section |5| 
The Kobulnicky fc Kewley [O III]/[0 II] diagnostic uti- 
lizes the instantaneous burst models shown in Figure^ 
Figure 13 shows the ionization parameter distribution for 
each sample assuming the instantaneous burst model. 
The logarithm of the mean ionization parameter is ~ 7.2 
for the NFGS and ~ 7.4 for our pairs sample, corre- 
sponding to a difference of ~ 8 x 10^ cm/s. Relative er- 
rors between ionization parameters estimated using the 
same calibration are typically << 0.1 dex. If contin- 
uous burst models are used, our ionization parameter 
estimates for the NFGS and the pairs are lowered by 
~ 0.3 dex, but the difference between the two samples 
remains ~ 8 x 10^ cm/s. We emphasize that it is not 
possible to determine whether differences in ionization 
parameter, burst scenario, or a combination of the two is 
responsible for the [O III]/[0 II] difference between the 
pairs and the fleld galaxies in Figure 0] 

The recent star formati on history of the galaxy pairs 
sample was investigated bv lBarton Gillespie et alJ l)2003l 
; hereafter BGK03). They showed that the equivalent 
width (EW) of Ha and other emission-lines correlate 
strongly with the inverse of the pair spatial separation 
(AD) or velocity separation. BGK03 used EW(Ha) 
with stellar population synthesis models to estimate the 
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Fig. 5. — The ionization parameter distribution for the galaxy 
pairs (soUd Hne) and the NFGS (dashed Hne), assuming an instan- 
taneous burst model. The logarithm of the mean ionization param- 
eter for each distribution is shown at the top. The logarithm of the 
mean ionization parameter is ~ 7.2 for the NFGS and ~ 7.4 for 
the pairs sample, corresponding to a difference of ~ 8 X 10® cm/s. 



time since the most recent burst began, or the "burst 
timescale". By matching the burst timescale with the 
dynamical timescale, BGK03 obtain a merger scenario in 
which the first close pass initiates a starburst, increasing 
EW(Ha). BGK03 note that this scenario is supported by 
the pairs data as long as the triggered bursts last longer 
than ~ 10* years. This timescale for star formation is 
consistent with our continuous star formation model in 
Figure 21 suggesting that the star formation history may 
contribute to the increased [O III]/[0 II] ratio seen in the 
galaxy pairs. 

Figure 01 shows that the metallicity predicted by the 
theoretical models is independent of burst scenario. The 
metallicity range spanned by the galaxy pairs is larger 
than the field galaxies; the pair sample contains more 
galaxies with relatively low metallicities (log(0/H)-|-12< 
8.8). Note that the [N n]/[0 II] ratio becomes insensitive 
to metallicity at low [N II] /[O II] ratios. In the following 
section we derive the metallicities for each sample and 
we investigate the cause of the metallicity distribution in 
Figure^ 

5. METALLICITY AND STAR FORMATION IN PAIRS 

We calculate metallicities for the pairs and the field 
galaxies using the strong optical emission-line ratios. 
There is a well-known discrepancy between metallicities 
calculated using theoretical strong-line diagnostics de- 
rivcd from photoionization models ('e.g..l McGaughl[l991t 
[Zaritskv et al. 1994; Kcwlcv ct al, 2002, and many oth- 
ers), diagnostics base d on electron-temperature mea- 
surements fe.g.. [ Kewlev fc EllisonI 120051 IGarnett et al.l 
I2nn4blal: iBresohn et al.ll2'ooir and diagnostics based on 
metallicities derived from a combination of the two meth- 
ods (e.g. Pcttini & Pagel 2004; Dcnicolo et al. 2002). 
This discrepancy manifests itself as a systematic offset 
with a maximum difference of ^ 0.4 dex in log(0 /H) + 12 
for star forming galaxies. The cause of this offset is un- 



known Csee IGarnett et al.ll2004bllat lStasihskall2005l for a 

discussion). Until this discrepancy is resolved, absolute 
metallicities derived using any method should be treated 
with caution. Fortunately, the difference between the 
various diagnostics is systematic. Therefore the error 
in relative metallicities derived using the same method, 
regardless of which me thod is used, is << 0.1 dex 
ijKewlev fc Ellisonll2005() . We therefore use metallicities 
derived using the same method only for relative compar- 
isons between samples or among galaxies within a sam- 
ple. 

We use a combination of the IKewlev et alJ lj2002f) 
"recomme nded" meth o d (he reafter KD02) and the 
[iTobulnick v Sz KewlevI l|2004j) recalibration of the 
IKewlev et all ()2002() [N II]/Ha-metallicity and R23- 
metallicity relations. These methods are based on a 
self-consistent combination of detailed current stellar 
population synthesis and photoionization models. These 
models su ccessfully reproduc e the H II region abundance 
sequence fDooita et al."2000') , the abundance sequence of 
star-forming galaxies (Kewle^^et al. 200 1a), and i ndivid- 
ual star clusters within galaxies ijGalzetti et al.l 12004). 
Our final metallicities are determined preferentially 
according to the following criteria: 

1. If [N II] A6584, [O II] A3727, Ha, and H/3 are 
available, and log([N II]/[0 II])> -1.0 we use 
the Balmer decrement to correct the [N II] /[O II] 
ratio for reddening and we apply the KD02 
[N II] /[O II] diagnostic to determine metallicities. 
The [N II] /[O II] ratio is extremely sensitive to 
metallicity but is insensitive to ionization parame- 
ter variations for log([N II]/[0 II])> -1.0. 

2. If [O II] A3727, [O III] A5007, Ha, and H^ are 
available, and log([N I I] /[O I I] )< -1.0, we use the 
Kobulnickv & Kewlev* l|2004f) recalibration of the 



Kewlev et al.l (2002) R23 -metallicity relation. 



3. If [N II] A6584, and Ha are available but 
[O III] A5007 and/or H/3 are not available we use 
the average of the Kobu lnickv fc Kcwlcv ( 2 00^ re- 
cafibration of the Ke wlev et al (2002) [ N II] /Ha- 
metall icity relation and the ICharlot fc Longhett) 
l)2001|) Case A [N II] /[S II] calibration. For high 
metalhcities ([N II]/Ha> 0.4) where [N II]/Ha 
becomes insensitive to metallicity we use the 
[N II] /[S II] calibration alone. We assume that 
the ionization parameter for these galaxies is the 
same as the sample mean ionization parameter. 
The mean ionization parameter is determined iter- 
atively using the reddening-corrected [O III]/[0 II] 
ratio as described in KD02. Deviation from the 
mean ionization parameter value is the domi- 
nant error source for metallicities calculated using 
[N II] /Ha and [N II] /[S II]. 

A total of 101/105 H II region-like galaxies in the 
NFGS and 153/153 pair members satisfy these criteria. 
We investigate the effect of residual AGN contribution 
on the derived metallicities in the Appendix. 

5.1. Luminosity-Metallicity Relation for Galaxy Pairs 

In Figure |^ we show the metallicity distribution for 
the pair members with projected separations 4 < s < 
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Fig. 6. — (a) The normalized metallicity distribution for the pair 
members with projected separations 4 < s < 20 kpc/h (red dotted 
Hne), 20 < s < 40 kpc/h (green soUd line), 40 < s kpc/h (blue 
dashed line), (b) the projected separation versus the metallicity for 
the pair members. The lowest metallicities occur predominantly in 
the close pairs. 



20 kpc/h (red dotted line), 20 < s < 40 kpc/h (green 
sohd hne), 40 < s kpc/h (blue dashed line). Figure 
plots the projected separation versus the metallicity. 
Galaxies with projected separations < 20 kpc/h (close 
pairs) have a smaller mean metallicity (8.85±0.02^) than 
galaxies with larger projected separations (8.97 ± 0.02). 

The metallicity distribution in Figure |^ for the close 
pairs is asymmetric and skewed towards lower metallici- 
ties compared to the more widely separated pairs. This 
difference is reflected in the luminosity-metallicity rela- 
tion. In Figure |7| we plot the luminosity-metallicity re- 
lation for the field galaxies (unfilled circles) and the pair 
members (filled circles) for the three ranges of projected 
separation. Panel (a) shows the least-squares best-fit to 
the L-Z relation for the field galaxies (dot-dashed line), 
and to the pairs in three groups of projected separation: 
4 < s < 20 kpc/h (red dotted hne), 20 < s < 40 kpc/h 
(green solid line), 40 < s kpc/h (blue dashed line). Pan- 
els (b-d) compare the L-Z relation for the pairs for each 
of the three pair groups with the L-Z relation for the 
field galaxy sample. Excluding the one outlier, the pairs 
with s > 40 kpc/h and with 20 < s < 40 kpc/h have 
similar L-Z relations to the field galaxies. This result 
is consistent with Marqugz et^ al. (2002) who find that 
the nuclear [N II] /Ha ratios for 13 non-disrupted galaxy 

Error is the standard error of the mean 
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Fig. 7. — (a) The least squares best- fit to the B-band lu- 
minosity metallicity relations for the field galaxies (dot-dashed 
line), and the pair members with projected separations 4 < s < 
20 kpc/h (red dotted line), 20 < s < 40 kpc/h (green solid line), 
40 < s kpc/h (blue dashed line), (b-c) comparison between the 
luminosity-metallicity relation of the field galaxies (unfilled circles) 
and the galaxy pairs (filled circles) for the three groups of projected 
separation in (a). The close pairs have a luminosity-metallicity re- 
lation shifted towards lower metallicities than the field galaxies or 
than the more widely separated pairs. 



pairs are similar to those of isolated galaxies; many of the 
Marquez et al. non-disrupted pairs are widely separated. 

Figure [Tji shows a striking difference; the close pairs 
have an L-Z relation shifted toward lower metallicities. 
For galaxies at the same luminosity, close pairs have 
lower metallicities than the field galaxies. Even more 
striking, there are 10 close pairs with metallicities lower 
than any field galaxy at the same luminosity. R-band 
luminosities produce identical results. 

This L-Z shift docs not appear to result from a dif- 
ference in luminosity distribution between the three pro- 
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Fig. 8. — The luminosity distribution for the galaxy pair mem- 
bers with projected separations 4 < s < 20 kpc/h (dotted line), 
and s > 20 kpc/h (solid line). The Kolmogorov-Smirnov two- 
sample test indicates that the two samples are consistent with a 
single distribution at the 72% confidence level. 



jected separation groups. We compare the luminosity 
distribution for the whole star-forming pairs sample (i.e. 
without luminosity cuts) in Figure |S1 . The close pairs 
(dotted line) do not lie at significantly higher luminosi- 
ties than the more widely separated pairs (solid line). 
The Kolmogorov-Smirnov two-sample test indicates that 
the two samples are consistent with a single distribu- 
tion function at the 72% confidence level. BGKOO calcu- 
lated the central burst strength su{t) of the galaxy pairs 
by fitting the observed B-R and Ha colors with stellar 
population synthesis models. Burst strength provides a 
measure of the current fraction of the central R-band 
flux originating from a new burst of star formation; pair 
members with recent bursts will have large central burst 
strengths. A total of 59 galaxies in our pairs sample have 
meas ured central burst strengths ijBarton Gillespie etHTI 
120031 hereafter BGK03). The recent central bursts con- 
tribute ^ 0.1 dex in R-band on average, with a range 
of 0.004 — 0.3 dex. Because the central burst strength 
is a measure of the central R-band flux originating from 
recent star formation, we can subtract the recent central 
star formation from the total R-band luminosity with 
knowledge of the slit-to-total R-band flux ratio and ask 
whether we can recover a "normal" L-Z relation char- 
acteristic of the underlying galaxy. Figure El shows that 
the corrected total R-band L-Z relation has the same L-Z 
shift as in Figure [7| 

Figure El is based on the assumption that the recently 
triggered star formation is located within the nucleus. 
However, even if we were to assume that the central burst 
strength is equal to the recent burst strength across the 
entire galaxy measured by the global R-band light, the 
R-band L-Z relations for the pairs and the field galaxies 
still do not agree. This exercise shows that the close pairs 
L-Z relation is not a result of a recent burst superimposed 
on a low luminosity, low mctallicity galaxy. 

Figures [3 and El show that the upper metallicity bound 
for the close pairs is lower than for the more widely sep- 
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Fig. 9. — (a) The R-band luminosity-metallicity relation of the 
field galaxies (unfilled circles) and the galaxy pairs (filled circles) 
for the three groups of projected separation as shown. The R- 
band luminosity for the pairs has been cor rected for the strength 
of the recent central burst of star formation lIBarton Gillespie et alj 
12003) . The dotted line shows the detection limit of the pairs sam- 
ple imposed by the recessional velocity selection criterion {cz > 
2300 km/s). While we must compare metallicities of galaxies at 
the same luminosity, we show the position of the field galaxies with 
cz < 2300 km/s for comparison. The close pairs have a luminosity- 
metallicity relation that is shifted towards lower metallicities than 
the field galaxies or than the more widely separated pairs, even 
after correction for the recent bursts of central star formation. 



arated pairs or the field galaxies. The wide pairs sub- 
sample (40 < s kpc/h) has 5/34 (15%) galaxies with 
metallicities log(0/H) -|- 12> 9.1. This fraction is com- 
parable to the fraction (4/21; 19%) of intermediate pairs 
(20 < s < 40 kpc/h) and field galaxies (8/43; 19%) with 
metallicities log(0/H) -I- 12> 9.1. In contrast, the close 
pairs sample has only 1/37 (3%) galaxies with metal- 
licities log(0/H) + 12> 9.1. Statistically, we would 
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expect 6-7 galaxies in the close pairs sample to have 
log(0/H)-|-12> 9.1 if there is no shift in the upper metal- 
licity bound with smaller projected separations. The 
deficit of close pairs with log(0/H) + 12> 9.1 as well as 
the shift in luminosity-metallicity relations for the close 
pairs leads us to conclude that the lower metallicities ob- 
served in the galaxy pairs probably result from a change 
in the gas-phase metal fraction in the nuclei of at least 
some of the close pair members. 

As a final check, we investigate the possibility of ex- 
tended triggered star formation in the galaxy pairs. Ex- 
ten ded star formation occurs in both s tarb urst galax- 
ies llRieke et al.l 1)19851) : ISchweize i* f 19861); see iSchweizeil 
l)2005|) for a revie w) and post -starburs t systeni s such 
as E+A g alaxies fFranx 1993; Caldw ell et all 119961: 
INorton et a l. 2001; Pracy et al. 2005). Our spectra rep- 
resent a luminosity- weighted mean of the individual star 
forming regions observed within the spectral aperture. If 
star formation is systematically triggered at increasingly 
large radii as the interaction progresses, the extended re- 
gions may dominate the light in the spectroscopic aper- 
ture at the late merger stages. In this case, our measured 
metallicities for the close pairs could represent the metal- 
licity in the extended regions rather than the nuclear star 
forming regions. If this process is important, we would 
expect to see an anti- correlation between aperture cov- 
ering area and metallicity, particularly for the close pairs, 
in the sense that the lowest metallicities would be mea- 
sured through the largest apertures. Our pairs spectra 
were obtained through a fixed-width aperture of 3 arcsec- 
onds. At the redshifts of our pairs sample, this aperture 
corresponds to a physical width of 0.01 - 2 kpc. The 
length of the slit varies from galaxy to galaxy, and at the 
distance of our sample, corresponds to a length of 0.01 - 
9 kpc, or 0.02 - 6.5 disk scale lengths. 

Figure E| shows the relationship between aperture 
length in units of disk scale length at the distance of 
each galaxy as a function of their residual metallicities 
from the field galaxy LZ relation (dashed line in Fig- 
ureO. We do not see an anti-correlation between metal- 
licity and aperture length. On the contrary, the lowest 
metallicities are measured through small-to-average sized 
aperture lengths (0.4 - 2.7 disk scale lengths, or 1-4 kpc). 
Similar results are found for aperture width and total 
aperture size. We conclude that the low metallicities in 
our close pairs probably do not result from star forma- 
tion being systematically triggered at increasingly large 
radii as the interaction progresses. 

We investigate the L-Z shift in terms of morphologi- 
cal and spectroscopic properties in Appendix None 
of these properties can explain the L-Z shift towards 
lower metallicities. We tentatively conclude that the 
lower nuclear metallicities observed in the close galaxy 
pairs result from a change in the gas-phase metal frac- 
tion. Because the lower-than-field metallicities occur 
predominantly in the close pairs, these lower metal- 
licites may result from an interaction-induced process. 
Interaction-induced gas infall from the outer, less en- 
riched regions of close pairs initiated by the first close 
pass is one possibility. Merger models predict that sub- 
stantial gas inflows occur towards the nucleus of disk 
^laxies undergoing interactions ("Barnes & Hcrnauistf 
| 199&: .Mihos fc Hernauist 1996: .B e kki fc Shiova 20Qlf 
ITono et a,lJl2004t ISnringel et a,ljl200l . Tn Sections 1^ 
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Fig. 10. — Nuclear aperture length (units of disk scale length) 
versus the residuals from the field galaxy LZ-relation (dot-dashed 
line in Figure for the galaxy pairs sample. The low- metallicity 
residuals lie to the right of the red-dashed line. The black line 
shows the position of the field galaxy LZ relation. The galaxies 
with large residuals do not have the largest aperture lengths. This 
result indicates that systematic triggering of starbursts at larger 
radii does not play a major role in the shift in LZ relation for close 
pairs. 



and l5.3l we investigate potential indicators of gas flows in 
our pairs sample. 

5.2. Tidally Triggered Central Bursts 

Theory predicts that gas flows in galaxy in- 
teractions trigger central bursts of star formation 
(jBarnes fc Hernauistill99^ . This infall is predicted to 
begin during the first close pass. As the galaxies move 
apart, EW(IIq:) decreases as the new population raises 
the continuum around Ha. Our sample may therefore 
select against galaxies that are enriched from the recent 
burst of star formation at large projected separations. 

We plot the nuclear metallicity versus the central burst 
strength in the lower panel of Figure^] The dotted line 
shows the least squares line of best-fit to the data. The 
metallicity and burst strength are anti-correlated. The 
Spearman Rank correlation test gives a correlation co- 
efficient of -0.43. The two-sided probability of finding a 
value of -0.43 by chance is 0.1%. If gas inflows have trig- 
gered the nuclear bursts, then those gas inflows may have 
carried more pristine gas from the outskirts of the galaxy 
into the nuclear regions, diluting the central metallicity. 
BGK03 showed that strong central bursts tend to oc- 
cur at close projected separations. Gas infall is expected 
to occur a fter two galaxies undergo their the first close 
encounter l|Barnesi r2002: Ba rnes fc Hernnuistl ir996:). 

The top panel of Figurellllshows that the pair member 
residuals from the mean field galaxy LZ relation (dot- 
dashed line in Figure Cj) are correlated with the central 
burst strength. Residuals are calculated as the short- 
est distance to the mean field galaxy LZ relation; posi- 
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Fig. 11. — Bottom Panel: Nuclear Metallicity versus central 
burst strength (sfl{t)) for the galaxy pairs. The dotted line shows 
the least squares line of best-fit to the data. The metallicity and 
central burst strength are anti-correlated. Top Panel: Residuals 
from the field galaxy LZ relation versus central burst strength. The 
strong correlation between LZ residuals and central burst strength 
indicates that the anti-correlation between metallicity and central 
burst strength is stronger than the strength that would result from 
a simple luminosity-star-formation and luminosity-metallicity scal- 
ing relationship. 
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tive residuals correspond to negative metallicities. The 
Spearman Rank correlation test gives a correlation coeffi- 
cient of 0.52. The two-sided probability of finding a value 
of 0.52 by chance is 0.007%. This result indicates that 
the anti-correlation between burst strength and metal- 
licity is over and above what one might expect from a 
normal luminosity-metallicity scaling relationship. 

In Figure 1121 we show the luminosity-metallicity re- 
lation for the pairs with central burst strengths mea- 
sured by BGKOO (solid symbols). We show the field 
galaxies for comparison (unfilled circles). BGKOO identi- 
fied six galaxies in the pairs sample that have high cen- 
tral burst strengths indicative of strong tidally triggered 
bursts (sfl(t) > 0.8). Of these six galaxies, five have 
luminosities between —21.5 < Mb < —18 and are dom- 
inated by star formation (diamonds in Figure IT^ . All 
five galaxies have a companion at close projected separa- 
tion s < 20 kpc/h. Their metallicities are all less than or 
equal to the lowest metallicity field galaxies at the same 
luminosity. The one remaining galaxy with a high central 
burst strength that does not satisfy our selection criteria 
is classified as 'Ambiguous'. This galaxy may be a com- 
posite so we have excluded it from our study. We note 
that this object is also in a close pair with s < 20 kpc/h, 
has a metallicity of log(0/H) + 12= 8.63 and luminosity 
Mb = —19.40, placing it firmly within the low metallicity 
region occupied by the five galaxies shown in Figure ISTTl 
The fact that all of the high central burst strength galax- 
ies have both close companions and predominantly lower 
metallicities than the field galaxies provides strong cir- 
cumstantial evidence that gas infall has played an impor- 
tant role in these low central metallicity galaxies. 



Fig. 12. — (a) The luminosity-metallicity relation for pair mem- 
bers with measured central burst strengths and projected separa- 
tions between (a) 4 < s < 20 kpc/h, (b) 20 < s < 40 kpc/h, and 
(c) 40 < s kpc/h. The five galaxies with strong central bursts (di- 
amonds; sji{t) > 0.8) are all in close pairs with metallicities lower 
or at the lower bound of the field galaxies. This result provides 
circumstantial evidence that gas flows have been important in the 
low metallicity close pairs. 



5.3. Blue Bulges 

Recent observations suggest that some bulges form 
within pre-existing disk galaxies ijKannappan et all 
^04). These bulges within disks may be linked to disk 
gas inflow and central star formation triggered by ei- 
ther intern al processes or galaxy mergers and interactions 
(jTissera et al. 2002). Active growth of bulges within 
disks can be identifled using B — R color profiles; galaxies 
which are actively growing bulges are likely to have bluer 
colors within their h alf-light radii compared t o the colors 
of their outer disks ijKannappan et all [2004 ; BGK03). 
We refer to bulges selected using B — R colors as "blue 
bulges" . 

BKG03 showed that the B ~ R color profiles of the 
six galaxies with strong central bursts discussed in Sec- 
tion |^1 also have blue central dips. Here we investigate 
whether the presence of blue bulges is correlated with 
metallicity. The B — R colors for the inner and outer 
disks are available for 87 of our galaxy pair members. 

In Figure 1131 lower panel we plot metallicity versus 
A{B — R) where A{B — R) is the difference between the 
outer disk (75%-light radius) color and the color inside 
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Fig. 13. — Lower Panel: Nuclear Metallicity versus A{B — R) 
for the galaxy pairs where A{B — R) equals the outer disk (75%- 
light radius) {B — R) color minus the {B — R color in the central 
half-light radius. Galaxies lying above the dotted line are likely to 
have blue bulges iKannaooan ct al. 20C3)- The dashed line is the 
least-squares line of best fit to the sample. Upper panel: Residuals 
from the field galaxy LZ relation versus A{B — R) for the galaxy 
pairs. The LZ residuals are weakly correlated with A{B — R). 
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the half-light radius (rg). The different symbols corre- 
spond to the three pairs groups with projected separation 
4 < s < 20 kpc/h (circles), 20 < s < 40 kpc/h (aster- 
ixes), 40 < s kpc/h (triangles). MetaUicity and A{B — R) 
are weakly anti-correlated; the Spearman Rank correla- 
tion coefficient is -0.29 and the probability of obtaining 
this value by chance is 0.7%. This correlation is domi- 
nated by the close pairs with low metallicities. 

Kannappan et al. classify galaxies with blue bulges 
as those galaxies with A{B — R) > (dotted line in 
Figure El • With this definition we see that galaxies 
with blue bulges cover the full range of metallicities in 
the pairs sample. However, galaxies with low metal- 
licities are a more likely to have blue bulges compared 
to higher metallicity galaxies; for metallicities below 
log(0/H)-|-12< 8.8, 33% of galaxies have blue bulges; for 
metallicities log(0/H) -I- 12> 8.8, 19% have blue bulges. 
In the blue bulge sample {A{B—R) > 0), a larger fraction 
of galaxies are in close pairs compared to the remainder 
of the sample (53% c.f. 34%). 

The top panel of Figure El shows that there is a weak 
correlation between the field galaxy LZ residuals and 
A{B — R). The Spearman Rank coefficient is 0.22 with 
a probability of obtaining this value by chance of 4%. 

In Figure 1141 we show the luminosity- metallicity rela- 
tion for the three ranges in projected separation. Galax- 
ies with blue bulges are shown with a large square out- 
line. Although most of the close pairs with lower metal- 
licities than the field galaxies have blue bulges, galaxies 
with blue bulges occupy the full range of luminosity, pro- 
jected separation, and metallicity in the pairs sample. 



Fig. 14. — (a) The luminosity-mctallicity relation for pair mem- 
bers with measured A(B — R) and projected separations between 
(a) 4 < s < 20 kpc/h, (b) 20 < s < 40 kpc/h, and (c) 
40 < s kpc/h. Pair members with blue bulges {A{B — R) > 0) 
are shown as squares. Although most of the close pairs with lower 
metallicities than the field galaxies have blue bulges, galaxies with 
blue bulges occupy the full range of projected separation and metal- 
licity in the pairs sample. 



6. DISCUSSION 

Our close pairs have lower metallicities, higher cen- 
tral bursts strengths and are more likely to be have blue 
bulges than galaxy pairs at wider separations. These 
low metallicity pairs are probably at a later stage in the 
merger process than pairs with higher metallicities. 

Smoothed particle hydrody namic merger simulations 
of 'B arnes &: HernauistI 1)199 6^ and iMihos &: HernauTs^ 
(1996) predict that a galaxy interaction inffuences the 
gas flow within a galaxy in two ways: (1) the tidal 
forces directly affect the gas by inducing gravitational 
torques which remove angular momentum from the gas, 
and (2) the gravitational influence of the companion de- 
forms the distribution of the stars and gas within the 
galaxy. The distorted stellar distribution changes the 
gravitational potential of the galaxy, which then affects 
the motion of the gas. The gas rapidly accumulates into 
the stellar arms that subsequently form shocked, dense 
filaments. In a dissipative medium, these filaments allow 
substantial streaming of gas towards the nucleus. The 
detailed theoretical response of t he gas in merging galax- 
ies was recently investigated by llono et all l)2004|) . lono 
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et al. predict that the average gas inflow rate to the 
central 1 — 2 kpc is 7 M0/yr with a peak infall rate 
of ~ 17 M0/yr. These simulations predict that the gas 
flows towards the central regions within 100 Myr after 
the initial collision, but before the two disks merge. Evi- 
dence for ionized gas flows in interacting galaxies has re- 
cently been found in Ha velocity maps of two galaxy pairs 
(jRampazzo et al.ll2005j) . and in radial ve locity cu rves in 
the interacting system Arp 194 by Marzia ni et al.l |2003). 
Cold gas flows have been observed in a number of merg- 
ing galaxies, inclu ding a galaxy pair separated by 29 kpc 
lichen et al.ll200l . 

BGK03 showed that the equivalent width (EW) of Ha 
and other emission-lines correlate strongly with the in- 
verse of the pair spatial separation or velocity separa- 
tion. By matching the burst timescale with the dynam- 
ical timescale, BGK03 obtain an interaction scenario in 
which the first close pass drives gas towards the central 
regions. This gas infall initiates a central starburst which 
increases EW(Hq!). As the burst ages and the projected 
separation increases, EW(Hq;) decreases as the new pop- 
ulation raises the continuum around Ha and enriches the 
gas. BGK03 note that this scenario is supported by the 
pairs data as long as the triggered bursts last longer than 
~ 10^ years. 

If a galaxy pair member is a late-type spiral with 
a strong metallicity gradient, then gas infall initiated 
by the first close pass could move less enriched gas 
from the outskirts of the galaxy into the central re- 
gions covered by our aperture. Massive late-type spi- 
rals like MlOl have large metallicity gradients; the H II 
region metallicities decrease by an order of magnitude 
from the inner to the outer disk (see e.g., IShield^lToM 
for a review). Metallicity gradients may result from 
one or a combination of: (1) radial variation of stel- 
lar yields caused by IMF differences between the spiral 
arms and the interarm regions (e.g. iGuesten fc Mezgeil 
Il982t): (2) a star forniation rate dependence on radius 
ijPhillipps &: Edmund£llT99l|) : (3) radial infall of primor- 
dial gas during disk formation (Matteucci fc FrancoisI 
[1989; PaEcl 1989; Edmunds & Grccnhow 1995-). Current 
chemical evolution m odels include a combin ation of these 
three processes ('e.g.. lChurches et a,1.ll2(¥)l . 

Not all galaxies have strong metallicity gradients, and 
the lack of a gra dient is oft en cited as evidence for ra- 
dial gas flows (see [Henry & Worthev 1999, for a review). 
Many barred spiral galaxies have weak er metallicity gra- 
dients than spirals of similar type (e.g.. lPagel et al.ll979l: 
iRov fc Walshll99'^ , suggesting that radial gas flo ws sup- 
press or change metallicity gradients (e.g., Ro berts et al.l 
[1979; Martin & Rov 1994; Rov & Walsh 199"^ (but c.f. 
IConsidere et al.ll2000r) . 

We can compare the fraction of infalling gas required 
to produce a shift in our metallicities with those pre- 
dicted by hydrodynamic merger models. The nuclear 
metallicity of our close pairs is ~ 0.2 dex lower than the 
metallicity of pairs at larger separations or fleld galax- 
ies at the same luminosity. At the mean redshift of our 
sample, our nuclear spectral aperture corresponds to the 
central ^ 1.5 kpc. Late-type spiral galaxies typically 
have metallicity gradients of ~ 1 dex from their nuclear 
to ou ter (8-9 kpc) H II regions (see IHenrv fc WorthevI 
119991 for a review). A 0.2 dex reduction in the nuclear 
metallicity in our pairs sample requires the gas within the 



central 1.5 kpc to be diluted by 50-60% with gas from 
the outer disk regio ns. Assuming a large initial central 
gas m ass of 10^ Mrr^ l) Sakamoto et al]ll999t iReean et alJ 
12001(1 . and an infall rate of 7 Mq/jv, to dilute the cen- 
tral gas mass fraction by 60% requires gas to infall over 
a timescale of ~ 9 x 10^ years. This infall timescale is 
consistent with the 1 x 10^ Myr timescale for gas infall 
predicted by the merger models. 

With a lower initial central gas mass of 1 x lO^M©, 
enough gas could accumulate from the outer regions to 
account for a 0.2 dex reduction in metallicity in only 
9 X 10^ years. The amount of dilution required is easily 
attainable in current merger simulations; merger mod- 
els predict that ~ 60% of the gas initially distributed 
throughout the progenitor disks may be driven into a 
small nuclear 'cloud' with dimensions of ~ 0.1 kpc dur- 
ing the merg ing process (Barnes & Hcrnauist 1991, 199^ 
lBarneslf2002^) . Evidence for such gas redistribution has 
been found by[Hibbard fc van Gorko m (1996). Hibbard 
et al. observed the distribution of cold and warm gas 
in galaxies along the Toomre sequence. They conclude 
that H I disks become more disrupted as the merger pro- 
gresses, from 60% of the H I mass being found in the 
main disks and bulges in each galaxy in the early merger 
stages to nearly 0% in the late merger stage. 

If the lower metallicities in our close pairs are a conse- 
quence of gas infall, these galaxies must be at a late stage 
in the infall process but at an early stage for chemical en- 
richment by triggered sta r formation . At 9. 4 x 10^ years, 
the merger simulations of llono et alJ l)2004D predict that 
the outer disk is no longer able to supply sufficient gas 
for infall, and the infall rate decreases. Observations of 
the metallicity gradients and the nuclear morphology in 
our close pairs could test whether the infall scenario de- 
scribed above is correct. 

We note that modeling gas flows in galaxy interactions 
is subject to uncertainties in our understanding of the 
kinematic behaviour of the gas. For example, a clumpy 
interstellar medium (ISM) may behave very differently 
in an interaction than the highly dissipative ISM that is 
traditionally assumed ( Charmandari s et al....200fl.'l . 

7. CONCLUSIONS 

We compare the nuclear metallicity and ionization 
properties of a large objectively selected sample of galaxy 
pairs and N-tuples with those of the objectively se- 
lected Nearby Field Galaxy Survey. We derive the first 
metallicity-luminosity relation for galaxy pairs and we 
compare this relation with the field galaxy luminosity- 
metallicy relation. We find that: 

1. Galaxy pairs with close projected separations s < 
20 kpc/h have systematically lower metallicities 
(~ 0.2 dex) than either the field galaxies or the 
widely separated pairs. There are 10 close pairs 
with metallicities lower than any field galaxy at 
the same luminosity. 

2. We find a tight anti-correlation between nuclear 
metallicity and central burst strength in the galaxy 
pairs; lower metallicity pairs have greater central 
burst strengths. This anti-correlation is stronger 
than that expected from the luminosity-metallicity 
relation alone. All five galaxies in the pairs sample 
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with extremely high central burst strengths have 
close companions and lower metallicities than the 
field galaxies. This result provides strong circum- 
stantial evidence that gas infall has plays an impor- 
tant role in the low metallicity galaxy pair mem- 
bers. 

3. We observe an anti-correlation between nuclear 
metallicity and A(i? — R). The luminosity- 
metallicity residuals indicate that this anti- 
correlation is stronger than the relation expected 
from the luminosity-metallicity relation alone. 
Galaxies which are actively growing bulges are 
likely to have bluer colors within their half-light 
radii compared to the colors of their outer disks. 
Pairs with low metallicities are more likely to have 
blue bulges compared to higher metallicity galax- 
ies; for metallicities below log(0/H)-Hl2< 8.8, 33% 
of galaxies have blue bulges while for metallicities 
log(0/H) + 12> 8.8, 19% have blue bulges. In the 
blue bulge sample {A{B — R) > 0), a larger frac- 
tion of galaxies are in close pairs compared to the 
remainder of the sample (53% c.f. 34%). This re- 
sult provides additional evidence for gas flows in 
the low metallicity close pairs. 

We use the [O III]/[0 II] vs [N II] /[O II] diagram to sep- 
arate ionization parameter from metallicity for our two 
samples and we show that the galaxy pairs have larger 
[O III]/[0 II] ratios on average than the field galaxies. 
This difference could result from either or a combina- 



tion of (1) a larger ionization parameter from a larger 
fraction of young hot stars in the central regions of the 
galaxy pairs, or (2) the galaxy pairs experiencing a longer 
duration burst of star formation than the field galaxies. 

Finally, we show that the combination of metallicity 
gradients typical of late-type spiral galaxies and gas in- 
fall as predicted by merger simulations can reproduce the 
reduced nuclear metallicities that we observe in our close 
pairs. In this scenario, galaxy interactions cause gas flows 
towards the central regions. The less enriched gas from 
the outer regions dilutes the pre-existing nuclear gas to 
produce a lower metallicity than would be obtained prior 
to the interaction. These gas flows trigger central bursts 
of star formation, producing strong central bursts, and 
possibly aiding the formation of blue bulges. Observa- 
tions of metallicity gradients and nuclear morphology in 
close pairs would test this scenario. 
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APPENDIX 
AGN CONTAMINATION 

In SectionOlwe classify the NFGS and galaxy pairs using their optical emission-line ratios. The lKewlev et af] l^2001a^ 
classification scheme removes galaxies that have an AGN contribution greater than ~ 20% to their optical line ratios. 
The potential residual 20% from an AGN could alter the line ratios and thus the metallicities. However, the majority 
of our pairs sample (99%), including all of the close pair outliers have metallicities determined using the [N II] /[O II] 
(77%) ratio or a combination of the [N II] /Ha and [N II]/fS II] (22%) hue ratio. A n AGN strengthens the [N II] and 
[S II I fines llBaldwin et al.lll981t iVeilleux fc Osterbrocklll987t iKewlev et al]l2001a(l but has little effect on the [O II] 
line (iFerland & Osterbrock'1986';'Ho et al."1993a'bVHo'"2005V Similarly, an AGN contribution increases the [N II] line 
relative to Ha (Baldwin et al. 1981; Vcillcux & Osterbrock 1987; Groves et al. 2004). Because both the [N II]/[0 II] 
ratio and the [N II] /Ha ratio increase with metallicity IfDenicolo et alJl2002t IKewlev fc DoDitall20"o3 iPettini fc Pagell 
|2004|) . an AGN contribution will cause our metallicities to be overestimated. 

The remaining 1% of our sample have metallicities determined using the R23 ratio. We use the R23 ratio when 
log([N II]/[0 II])< —1.0. This limit corresponds to a metallicity of log(0/H) -I- 12< 8.4. In this low metallicity regime, 
the R2.'^ ratio also rises with increasing metallicity. An AGN contribution increases [O III] relative to [O II] or H/? 
ijBaldwin e t al. 1981: V eilleux fc Osterbrockll987tllIao et al.ll2005t 11101120051) . thus increasing the R23 ratio and causing 
the the resulting metallicity to be overestimated. 

Therefore, a systematic contribution from a low-level AGN to the close pairs sample would cause an artificial 
increase in our metallicities which would shift the luminosity-metallicity relation in the opposite direction to that 
which is observed. We conclude that an AGN contribution can not cause the observed shift in L-Z arelations toward 
lower metallicities for the close pairs. 

MORPHOLOGICAL AND SPECTROSCOPIC EFFECTS 

In this section we investigate whether real morpholog ical o r spectroscopic effects might account for the shift in the 
luminosity-metallicity relation for close pairs. In Figure lB15b . we plot the residuals of the galaxy pairs from the mean 
luminosity-metallicity relation of the field galaxies (dot-dashed line in Figure Ej) versus the mean equivalent width 
(EWm) of the combination of lines used to determine metallicities. For most galaxies in the sample, metallicities were 
determined using the [N II] /[O II] ratio with reddening correction using the Ha/H/3 ratio. In this case, EWm is the 
mean equivalent width of [N II], Ha, [O II], and H/3. Galaxies with low EWm have the most uncertain metallicities. 
The red dashed line marks the low metallicity outliers of the mean field galaxy LZ relation. The EWm for the outliers 
is high, indicating that the metallicities derived for them are not systematically affected by signal-to-noise issues. 
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Fig. B15. — The residuals of the galaxy pairs sample from the mean luminosity-metallicity relation for the field galaxies (Figure|3 versus 
(a) The mean equivalent width of the lines used to calculate metallicities for each galaxy, (b) de Vaucouleurs Morphological Type, (c) Bar 
Class. Galaxies with strong bars, weak bars, and no bars are classed as 2,1,0 respectively. We have shifted the close and widely separated 
pairs +0.1 and -0.1 from these bar classes respectively so that their distribution can be viewed more readily. The black line shows where the 
pair members should lie if they were distributed uniformly about the least-squares fit to the field galaxies luminosity-metallicity relation. 
The red dashed line shows the position of outliers in Figure 171 Positive residuals correspond to lower-than-field metallicities. 



IKannappan et af] lj2004|) analysed the morphologies of the pairs sample. In Figure IB 15b . we show the relationship 
between de Vaucouleurs morphological type and the residuals. The outliers do not have significantly different range 
of types than the remainder of the sample. 

Finally, we compare the presence of a bar w ith the residuals in Figure IB15b . The presence of bars in the pairs 
sample was studied by ijKannappan et aDl2004|) . Galaxies with a strong bars, weak bars, and no bars are classed as 
2,1,0 respectively. The close and widely separ ated p airs are shifted -1-0.1 and -0.1 from these classes so that their 
distribution can be viewed more readily. Figure IB 15b shows that the fraction of bars in the outliers is low. Thus, the 
presence of bars cannot account for the shift in the LZ relation for the close pairs. 
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